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ABSTRACT 

A diffusion-limited annihilation process, A+B — > 0, with species initially 
separated in space is investigated. A heuristic argument suggests the 
form of the reaction rate in dimensions less or equal to the upper critical 
dimension d c = 2. Using this reaction rate we find that the width of the 
reaction front grows as t 1 / 4 in one dimension and as t 1 / 6 (lnt) 1 / 3 in two 
dimensions. 

P. A. C. S. Numbers: 82.20.-w, 05.40.+j., 02.50.-r, 82.70.-y. 
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In this paper, we investigate the kinetics of diffusion-limited two-species annihilation 
process, A + B — > 0, in which collisions between two distinct species A and B lead to 
the formation of inert reaction products. Classically, such a process can be described by 
reaction-diffusion equations [1], 

da ^d 2 a db ^d 2 b , . 

— =D r — = D r (1) 

dt dx 2 ' dt dx 2 y } 

Here a = a(x,t) and b = b(x,t) are the local concentrations of A- and S-species, re- 
spectively, D is the diffusion coefficient which is assumed to be identical for both species, 
and r = r(x,t) the reaction rate. In the mean-field approximation this reaction rate is 
proportional to concentrations of reactants, r = kab. For sufficiently low dimensions, the 
diffusion mechanism is not efficient enough and fluctuations in the densities of diffusing 
reactants will result in a different form for the reaction rate and dimension-dependent 
kinetic behavior at long times [2] . 

Several recent studies have been focussed on the situation in which species are initially 
separated in space and hence they react in a confined region called the reaction front. This 
front plays an important role for a variety of physical and chemical processes [3]. At the 
mean-field level, a scaling theory has been developed [4] which predicts that the width 
of the reaction front, w, grows in time as t 1 / 6 in agreement with experiments [5] and 
simulations [5-7]. 

Since the upper critical dimension for this process appears to be d c = 2 [7] , a departure 
from the value a = g for the width exponent is expected for d < 2. Numerical simulations 
of one dimensional systems suggest the value a as 0.3 [7-9]. However, an exact value of the 
width exponent in one dimension is unknown. In two recent works it was argued that ol = \ 
in one dimension. Both these studies are based on analysis of more tractable systems: the 
first paper [10] examines the case in which one reactant is static and the other diffuses 
while the second paper [11] treats the behavior of the front in the steady state reached by 
imposing the antiparallel current densities of A and B species at x = +oo and x = — oo, 
respectively. The latter steady state version of the original model was introduced in [12] 
and further explored in [13]. 

In the present study, we investigate the original transient problem of initially separated 
equally mobile species. Our analysis is based on the rate equations approach but with a 
modified reaction rate. In one dimension, we find a = \ in agreement with previous studies. 



In two dimensions we derive a logarithmic correction to the mean-field behavior: w ~ 
t 1 /6( lnt )i/3_ while it may be difficult to confirm this logarithmic correction numerically, 
a significant departure from the mean-field behavior, w ~ t 1 / 6 , has been observed [7]. 

Let us first estimate the reaction constant k in the mean-field relation r = kab. It is 
reasonable to assume that the reaction constant depends on the diffusion constant D and 
the radius of particles R which, for simplicity, is assumed to be the same for both species. 
Simple dimensional analysis gives k ~ DR d ~ 2 . On physical grounds, the reaction constant 
must be an increasing function of the radius R and hence the previous expression for k is 
valid only for d > 2. To find the reaction constant for d < 2 we should assume that k does 
not depend on R but instead it depends on concentrations, k = k(D, a, b). A dimensional 
analysis gives k ~ D(a + 6)~ 1+2 / d and therefore 

r ~ Dab(a + b)- 1+2/d . (2) 

In fact, any symmetric function r = r(a, b) of densities a and b with degree of homogeneity 
1 + 2/d is equally possible. However, in finding scaling exponents the precise form of the 
reaction term is unimportant and for concreteness we choose the form given above. Note 
that this simple dimensional argument also correctly predicts the value of the upper critical 
dimension, d c = 2. 

This type of argument was first applied by Toussaint and Wilczek [14] to homogeneous 
single-species annihilation, A+A — > 0, and led to the qualitatively correct long-time kinetic 
behavior: a ~ (Dt)~ d / 2 for d < 2. In contrast, for homogeneous two-species annihilation, 
A + B — > 0, this argument predicts the incorrect value of the upper critical dimension 
(the correct value is d c = 4) and the incorrect kinetic behavior (t~ d l 2 instead of t _d / 4 , 
see [2,14]). Because of this difficulty there were no attempts to apply the argument to 
inhomogeneous two-species annihilation. It can be readily seen that a large-scale spatial 
organization of reactants [2,14] for homogeneous two-species annihilation leads to this 
difficulty. Indeed, the reaction takes place only near boundaries between different-species 
domains. The fraction of space occupied by these reaction zones tends to zero as t — > oo 
[15]. Hence modeling the reaction term by a function of averaged concentrations a(t) and 
b(t) is inappropriate since it gives the same reaction rate inside single-species domains 
where reaction does not proceed and in the boundaries. However, for inhomogeneous two- 
species annihilation, the behavior in reaction zone where the reaction does proceed is of 
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interest. Therefore, a description of reaction zone by using a (modified) reaction rate which 
is a function of local concentrations a(x, t) and b(x,t) seems reasonable. 

Consider now a one dimensional inhomogeneous system in which A particles are uni- 
formly distributed to the right of the origin and B particles are uniformly distributed to 
the left of the origin, both with equal concentration cq. The difference a(x,t) — b(x,t) 
satisfies the diffusion equation which, subject to given initial conditions, has the solution 



a(x, t) — b(x, t) = Co erf ^ — j== ^ . 



(3) 



In the long-time limit, the scaling form 

a(x,t) ^t-^A(xt- a ), b(x,t) ~ t" 7 i3(xt- Q ), r(x,t) ~t- p 1l{xt- a ), (4) 

is expected. Since r ~ (a + b)ab in one dimension (see Eq. (2)), the scaling form of 
Eq. (4) implies the scaling relation f3 = 37. Further, by inserting Eq. (4) into Eq. (2) and 
asymptotically balancing various terms one finds another scaling relation 7 + 2a = 37. 
Finally, assuming the reaction zone is short compared to the diffusion length \^Dt and 
making use of the exact solution (3) one gets a — b ~ cox/VnDt ~ t a_1 / 2 . This gives the 
last scaling relation a + 7 = |. Combining these scaling relations we find the exponents: 
a = 7 = |,/3=|. Note that the width indeed increases slower than the diffusion 
length, w ~ t 1 / 4 <^ t 1 / 2 , thus confirming our assumption. Rewriting the scaling ansatz in 
dimensionless form, we arrive at 



r 2 \ 5 / r 2 \ 4 / r 2 



a ( x >*)-{jji) ^(0, 6 (^*)^^J r{x,t) ~ D n(0, (5) 

with the scaling variable £ = x (c 2 / 'Dt) 1 4 '. 

Substituting Eqs. (5) into the governing equations (1) and making use of Eq. (3) we 
arrive at still intractable nonlinear ordinary differential equation. However, if we focus 
on the region far enough from the origin, £ ^> 1, some progress can be made. To find 
the scaling form of densities we have to choose an appropriate form of the reaction rate. 
Let us mention that the reaction rate given by Eq. (2) at d = 1, r ~ Dab(a + 6), has 
an additional advantage of being linear in concentration which is small compared to the 
other, e. g., r ~ Da 2 b oc b if b a. This physically appealing behavior suggests to use 
this particular form. 



So, let us examine Eqs. (1) with the reaction rate 



r = KDab(a + b), (6) 

where k is a dimensionless constant. In the region i^Dt/c^) 1 ^ C 1 <C \fT5i where the 
A-species dominates Eq. (3) simplifies to a ~ c$x / which can be rewritten as 
A{^) ~ £/ 1 ypH in the scaling form. Then the second Eq. (1) asymptotically reduces to 

d 2 B K 



n 



(7) 



By applying the WKB procedure [16] one finds the asymptotic solution to Eq. (7), 

1 



S(0 ~ ^ ex P 




(8) 



which is valid for £ 3> 1, i. e., for x 3> (Dt/cl) 1 ^ 4 . Combining Eqs. (5), (6), and (8) we 
arrive at the final asymptotic expression for the reaction rate: 



' 2 \ 9 / 8 



DtJ 



x ' exp 




k 2 / c o 

X ' 



2 \ 1/2' 
0_ 

Dt 



(9) 



for (Dt/cl) 1/A < x < v^Dt. 

Let us compare our prediction for the tail of the scaling distribution, Eq. (9), with 
other studies. In a recent paper [8], the validity of the scaling distribution has been 
challenged and a description based on a continuous spectrum of scales between t 1 / 4 and 
t 3 / 8 has been proposed. In particular, Araujo et al. have fitted the tail of the reaction 
front by the exponential form: r ~ exp [— |ic | / 1 3 / 8 ] . In a very recent study [17], Cornell has 
reconsidered the problem. Based on extensive simulations he has found that the dynamical 
scaling appears to hold, contrary to the claim of Ref. [8]. He has obtained a 0.28 for 
the width exponent and observed that the reaction rate profile has a Gaussian form. Thus 
our results well agree with the most extensive present simulations [17]. 

Let us now turn to the two dimensional version of the problem. The same line of 
reasoning gives the reaction rate r ~ Dab and hence the mean-field behavior could be 
expected. However, d = 2 is the marginal dimension, logarithmic corrections to mean-field 
expression r ~ Dab may arise, and therefore one has to be more careful. To assess the 
validity of the mean-field expression for the reaction rate let us start from more fundamental 



relation, r ~ b/T, where T is a time interval in which a B-particle typically collides with 
some A-particle. Consider for simplicity the region in which A species dominates. In 
a reference frame at rest with respect to an arbitrary "target" B-particle, the density 
a(x + r, t + t) of A-particles is governed by the diffusion equation 

da(x + r, t + t) 



dr 



2DA r a(x + r,t + r). (10) 



We should solve this equation subject to the initial condition a(r = 0) = a(x,t) (since 
the local density is a(x,t)) and the adsorbing boundary condition a(|r| = 2R) = for 
r > 0. When Dt 3> R 2 , an approximate solution to Eq. (10) with prescribed initial and 
boundary conditions can be readily found by using a quasistatic approximation (see, e. g., 
[18]). This solution reads 

, , \n(r/2R) 

^ + r,« + T) e! a(»,«) 21n ^ ) . (11) 

Note that the "local" density given by Eq. (11) is changed on the length scale of the 
order V Dt which has to be small compared to the length scale of the "external" density 
a(x,t), i. e., \TD~t <^ w. We will assume that both inequalities R 2 <^ V DT <^ w hold 
and therefore Eq. (11) may be applied to estimate the collision time T. From the final 
expressions for T and w it is straightforward to verify that previous inequalities are indeed 
satisfied in the long-time limit thus providing the check of self-consistency. 

The collision time T may be evaluated by computing the flux to the circle of radius 
2R and then by equating the flux to the unity: 

SwDR jsjM^M±ll dT = 1 . (12) 



By combining Eq. (12) and Eq. (11) one gets the final estimate, T ~ In (l/2ixaR 2 ) /2nDa. 
Thus we obtain the following expression for the reaction rate, 



r = 2nDab 



1 1 

+ 



(13) 



_ln(l/27rai? 2 ) In (l/27ibR 2 ) _ 

The second term is added to keep r = r(a, b) symmetric. In the region x ^> w which we 
have considered the second term is small and thus never appeared. Similarly, in the region 
—x — w the first term is negligible. Since the reaction rate must be symmetric function 

6 



of the densities, Eq. (13) provides a proper combine approximation in both tail regions. It 
is therefore reasonable to assume that Eq. (13) is in fact accurate everywhere in reaction 
zone. So instead of the mean-field expression for the reaction rate in two dimensions, 
r = Dab, a refined argument gives a logarithmically corrected value (13). 
Now let us try a scaling ansatz 



<x,t) ~ ( |0 3 LM(0, b(x,t) ~ (||) 3 L"B(0, 
r(x,t)^D(-|) 3 L^(0, ^^(^'^ 



(14) 



where we have used a shorthand notation L for the logarithmic factor L = In (Dt/cg-R 6 ). 
Eq. (14) is just the standard mean- field scaling ansatz [4] modified by logarithmic cor- 
rections. A procedure identical to the one used in the one dimensional situation gives 
three scaling relations between "logarithmic" exponents A, fi, and v from which one finds: 

\ = V=-\,\L=\. 

The main lesson from this study is the following: Even for sufficiently low dimensions, 
d < d c , when the strong correlations destroy the mean-field approximation, it is still 
possible to use mean-field like diffusion-reaction equations with modified reaction rates. 
At the critical dimension, the present heuristic method gives subtle logarithmic corrections 
which seem hard to obtain by more solid theoretical approaches. 

Let us finally mention that the more general higher order annihilation reaction process, 
mA + nB — > 0, with initially separated species can be analogously investigated by the 
present approach. Since the critical dimension for this reaction process is d c = m+ ^_ 1 
(see, e. g., [11]), deviations from the mean- field predictions are expected only for the 
classic case (m, n) = (1, 1) [when d = 1 and 2] and for (m, n) = (2, 1) [when d = d c = 1]. 
Consider the latter case and assume that point-like particles undergo a random walk on 
the linear lattice with spacing Ax. [Even in one dimension we must keep the lattice 
spacing finite since otherwise the reaction with three-particle interaction is absent.] Focus 
again on the region where A-species dominates and estimate the collision time T. Note 
that the collision event, i. e., the situation when two A-particles will be simultaneously in 
the site with a target S-particle, may be described as usual two-body collision between 
an immobile target particle and imaginary particles corresponding to pairs of A-particles 
and diffusing on the two-dimensional square lattice with the same spacing Ax. A more 



detailed description of this correspondence is given in Ref. [19] which considers a related 
single-species aggregation process, A + A + A — > A. Then up to the replacement of the 
one-dimensional density a by the two-dimensional density of pairs a 2 and the spacing 
Ax by the radius R one can use previous results (see Eq. (10)-Eq. (12)) and find the 
collision time: T ~ In (l/2ira 2 Ax 2 ) /2ixDa 2 . So we obtain r ~ Da 2 b/ ln(l/aAx) while 
the mean-field approach would give r ~ Da 2 b. This results in logarithmic corrections 
to the mean-field behavior. Repeating the steps employed for the previous treatment of 
two-particle annihilation we arrive at the scaling form 

a(x, t) ~ (L/t) 1 / 4 ^), b(x, t) ~ (L/O^BCO, r(x, t) ~ r^L-^n{0. (15) 

Here ^ = a; (cgZ-Dtiv) 1 ^ 4 is the scaling variable and L = In (.Di/cgAx 4 ) the logarithmic 
factor. Thus we conclude that the width of the reaction front for the three-particle reaction 
process, 2 A + B — > 0, scales as t 1 / 4 (lnt) 1 / 4 in one dimension. 

I am thankful to S. Redner for discussions and a critical reading of the manuscript 
and to E. Ben-Nairn for correspondence. This work has been partially supported by grants 
from ARO and NSF. 
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